Effect of the intrinsic Josephson coupling on the pancake lattice in layered 

superconductors. 
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We study the pancake lattice system induced by application of transverse magnetic filed to the 
layered superconductor. A simple statistical field theory for the pancake lattice is derived. It 
incorporates effects of the magnetic interaction between pancakes as well as those of the interlayer 
Josephson coupling. The proposed description enables us to estimate the pancake crystal melting 
temperature from above. Also, it allows us to compare directly the effects of the magnetic and 
Josephson energies on the statistical properties of the pancake lattice. We demonstrate that even 
for such an anisotropic material as BSCCO the pancake interaction induced by the Josephson 
coupling is of the same order as the magnetic interaction. 



stances it is possible to drop the Josephson coupling from 
the description and when it cannot be omitted. 

The paper is organized as follows. In Sect. II start- 
ing from the vortex Hamiltonian for the anisotropic su- 
perconductor we reproduce the well-known result for the 
magnetic potential energy of the pancake lattice. In Sect. 
Ill the energy due to the Josephson coupling is derived. 
In Sect. IV we compare Josephson and magnetic interac- 
tion energies. The purpose of the latter comparison is to 
determine when it is permissible to ignore the Josephson 
mechanism and when it is not. In Sect. V we apply the 
derived Hamiltonian to the problem of pancake lattice 
melting. Sect. VI contains conclusions. 



I. INTRODUCTION 

Magnetic field penetrates a sample of layered supercon- 
ductor forming pancake vortices in every layer of the sam- 
ple. Calculation of statistical properties of the pancake 
matter is a very complicated question for several reasons 
of which the following is the most important. The pan- 
cakes interact through two mechanisms, magnetic and 
Josephson. While magnetic mechanism is reduced to 
pairwise interaction of the pancakes, potential energy due 
to the interlayer Josephson coupling cannot be brought 
down to such simple form. Available methods lead to 
a non-local description [1] or rely on the duality trans- 
formation [2, 3] which maps the pancake system onto 
a system of quasi-two-dimensional Coulomb-like plasma. 
Both approaches are technically complicated and difficult 
to generalize. 

Alternatively, many authors neglect Josephson mecha- 
nism assuming that for extremely anisotropic supercon- 
ductors its influence is small compared to the effects of 
the magnetic interaction [4, 5, 6, 7, 8]. 

In this paper we show that under certain circumstances 
it is possible to include Josephson coupling into consid- 
eration with the help of general physical devices. At suf- 
ficiently high field the Josephson and magnetic interlayer 
coupling get smaller than the intralayer interaction. In 
such a situation the system is viewed as a collection of 
weakly coupled 2D layers. If pancakes in every layer form 
a structure with short range crystalline order both mag- 
netic and Josephson mechanisms could be accounted for. 
The resultant model is a version of the sine-Gordon field 
theory. It could be analyzed with the help of usual field 
theoretical tools. 

The proposed approach allows one to find an upper 
bound for the pancake lattice melting temperature. An- 
other interesting application is that it becomes possible 
to compare contributions of the magnetic and Josephson 
interactions to the statistical properties of the pancake 
matter. Thus, we could determine under what circum- 



II. MAGNETIC COUPLING 

We start with the vortex Hamiltonian H [9, 10]: 
H^^ I ^/C,(q)|z;,(q)|'+/C||(q)|«||(q)|'}, (1) 
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The symbol ■ • ■ stands for J (2'K)^^(i?'q . . .; vector field 
V denotes vorticity. The component in the nth super- 
conducting layer coincides with the pancake density in 
that layer: 



v^{v, n) = p{r, n) S^'^'^ (r - r^(n)) 



(4) 



where r^(?T.) is the position of the fith vortex in the nth 
layer. 

The parallel components of the vorticity V|| correspond 
to the Josephson vortices between superconducting lay- 
ers. The Josephson vorticity between nth and (n — l)th 
layer is equal to: 

v|| (r, n) = ^ ^{r^(n) - r^(n - l)}Sf,{r), (5) 



2 



W= (6) 

/ _ _ 1) _ _ _ 1)]^-, _ 

Jo 

Here is the spacing between neighboring layers, r is 
2D coordinate vector. Function is defined in such 
a way that it is locahzed on a straight Une connecting 
points r(n) and r(n— 1). Physically, it is non-zero on the 
Josephson vortex which connects a pancake in nth layer 
and another pancake in (n — l)th layer. 

Since a vortex cannot terminate inside the supercon- 
ductor we have 

divv = 0, (7) 
or, equivalently, in Fourier space 

QzVz = (8) 

Let us rewrite the kernel JCz in the following manner: 
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Note, that /C 



mag 



contains no Ac- In other words, it does 
not depend on the Josephson coupling between layers. It 
describes magnetic interaction between pancakes. On the 
contrary, kernel 6K,j vanishes when the Josephson energy 
vanishes. 

Now, using 3D-transversality of the vorticity (eq.(8)), 
we establish the following: 
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This allows us to rewrite the Hamiltonian (1): 
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In other words, only the 2D-transverse part V||t of the 
Josephson vorticity V|| contributes to the interaction en- 
ergy. 

Magnetic part of the energy is the easiest to handle. 
Imagine that in every layer the pancakes form 2D crys- 
tal with short range hexagonal lattice characterized by 
the lattice constant uq. We assume that it is possible to 
choose the length scale ^ ^0 in such a way that within 
^0 the lattice translational invariance is undisturbed. Al- 
though, we do not need the exact value of for our 



calculations, yet, our ability to set such a scale is imper- 
ative. 

If the above prerequisite is met a pancake lattice in 
nth layer could be described locally by a 2D displacement 
vector field u„. This field does not change much on the 
scale ^0- The pancake density is expressed as: 

v^r, n)=poU + div2D u„ + ^ e*iiir+iqiiu„ j . (le) 

The sum is upon vectors of the reciprocal lattice. The 
2D divergence is defined as: 



dUa: , dUy 
dlV2DU = ^— + — ^. 

OX oy 
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The quantity po denotes the average density of the pan- 
cakes: 



B 



(18) 



Observe, that two first terms of (16) vary slowly on 
the scale of the lattice constant ao while the third term 
quickly oscillates as function of r: 



Vz = P + Sp, 
p = Po (1 + div2£) u„) , 
5p = poY^ e^iii'+'^ii^n. 



(19) 
(20) 
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We have for the magnetic interaction energy: 



Wmag- 



J2 f^m^g{^-r',n-m)v^v'^= {22} 
-ff- Z] J , ''^mag(r - r', n - m) {pp' + 5p5p') , 



where the symbol • • • stands for j cPr . . .. The term 
proportional to pp' describes the statics of the longitudi- 
nal displacements. We will not be interested in them for 
they are suppressed [11]. The remaining parts are: 



^mag = '^Pldl X (23) 



X] /C„,ag(r-r',n-m) ^ e*'iii('-'"')+*'»ii 



The notation 



(24) 



was adopted in the formula above. To avoid clutter we 

omit the subscript 'T' which denoted the transverse dis- 
placement mode: we write here u„ rather than u„t- We 
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will continue to do so, yet, one has to remember that the 
displacement vector field we consider is pure shear. 
The first term in the above expression gives usual Cqq: 



66 
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To proceed further with the second term we need to cal- 
culate /Cn 



-mag- 



ciqzdoU 
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If XabQ\\ ^ 1 which is true for B > Hci the above expres- 
sion for /Cmag could be simplified: 



mag ~ ^2 „2 \ "'c 



1 f _^ e-<l\\<ic\n\ 



Thus, we get for the second term in (23): 



(27) 
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^-q]]dc\n-m\ r 



„iq||(u„-u„) 



Here we neglected the weak dependence of u(r + Ar) on 
Ar since |Ar| ~ oq <C ^o- Finally: 
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Here qo is the magnitude of the smallest reciprocal lattice 
vectors: 



90 
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(30) 
(31) 



Properties of the Hamiltonian similar to (29) was inves- 
tigated in [4, 5, 6, 7, 8]. 

At finite temperature the above expression could be 
simplified. It is enough to notice that at T > contri- 
bution of a term cxp(?'q|| (u„ — u^)) is proportional to 
exp(— (7^(u^J/2). Thus, when fluctuations of the lattice 
are substantial it is permissible to retain in (29) lowest 



|q||| terms only. Consequently, it is not necessary to sum 
over all possible qy in the above formula. It is sufficient 
to keep only six terms corresponding to six elementary 
reciprocal lattice vectors whose absolute values are equal 
to qo- Namely, the last term of (29) might be written as: 
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-qodo\r. 



cos(q|i(u„ - Um))- 



Our model becomes a vector version of a sine-Gordon 
field theory. 

III. JOSEPHSON COUPLING 

Now we add the Josephson coupling to our Hamilto- 
nian. This means that we have to include the Josephson 
vorticity as well. At first it seems like an impossible task 
since we have to account for the Josephson vorticity fluc- 
tuations. We will demonstrate that under rather broad 
conditions these fluctuations could be included by simple 
renormalization of the Josephson coupling parameter. 




FIG. 1: A fluctuating Josephson vortex (solid line) connects 
two pancakes, one in the nth layer (closed circle), another 
in the (n — l)th layer (open circle). The deformed vortex 
is decomposed into a straight line and a loop (both broken 
lines). The characteristic size of the loop is fioop- 

On fig.l an example of Josephson vorticity fluctuation 
is shown. In the absence of fluctuations two pancake 
vortices in neighboring layers would be connected by a 
straight Josephson vortex. Fluctuations create an out- 
growth on the vortex line. The size of the outgrowth is 
Cioop • The deformed vortex is represented by a solid line. 
The deformed vortex could be "decomposed" into a vor- 
tex loop and a straight vortex connecting two pancakes. 
Thus, the problem of the Josephson vorticity fluctuations 
could be reduced to the problem of the thermally induced 
vortex loop gas. 
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Thermal fluctuations infuse superconductor with loops 
of Joscphson vortices of typical size ^loop- This size ^loop 
diverges at the superconducting transition temperature 
Tc and vanishes at T = 0. It is finite at Tc > T > 0. 
Loosely speaking, ^loop characterizes a spacial scale above 
which "mean-field" Josephson coupling could be defined 
and fluctuations are unimportant. Therefore, if we are in 
the regime 



6oop < ao 



(33) 



the fluctuations might be omitted. In such a case fluctu- 
ations result in small vibrations of the .Josephson vortex. 
The only consequence of the fluctuations would be renor- 
malization of eff'ective value of the Josephson coupling, 
which is equivalent to renormalization of the penetration 
depth from its bare value Af to experimentally measur- 
able value of Ac. 

Of course, it is desirable to generalize the argumenta- 
tion beyond (33) . How this could be done is discussed in 
Appendix A. 

Once we settle the issue of the fluctuations we could 
proceed with the derivation of the Josephson contribution 
to the energy of the pancake system. 

As in the previous Section we assume that in ev- 
ery layer the pancakes form a 2D structure with (at 
least) short range crystalline order. Remember also that 
the longitudinal displacements of the pancakes are sup- 
pressed. In absence of the longitudinal displacements the 
Josephson vorticity is a periodic function of vector s„ de- 
fined as (fig. 2): 

s„(r) = u„(r) - u„_i(r). (34) 

Two-dimensional vector field vy (r, n) could be expressed 




where the smooth part is equal to: 
V|| = ^g(s„). 

Vector function g is defined like so: 
g(s) = s - l(s). 



(36) 



(37) 



The lattice vector I(s) is chosen to deliver minimum to 
the expression |s — 1| (see fig. (3)). Function g(s) is invari- 
ant under lattice translations. For any s vector g always 
lies within the primitive lattice cell. Physically. eq.(36) 
means that when two lattices are shifted only slightly 
with respect to each other V|| = pas/dc- However, when 
s does not fit into the primitive cell the Josephson vor- 
tices rearrange themselves to minimize their own length. 
In this case |g| < |s|. For example, if s is the lattice 
vector then g = 0. 

Next, we calculate the oscillating part. Full Josephson 
vorticity field between nth and (n — l)th layers is given by 
eq.(5) with rn{n — l) =1,^-1- u„_i (I^ is the lattice vector 
which shows the undisturbed position of the /ith pancake) 
and r^(n) — r^(n — 1) = g. The Fourier transform of the 
vorticity vector field due to a single vortex n is: 



V||,^(q||,n) 



J-c Jo 



;(s„) 



The total vorticity is: 



^giq||(u„_i 



-i)e*i|i'^. (38) 



= X!'^llA'(1ll'")• 



(39) 




FIG. 2: The Joscphson vorticity generated by the relative 
displacement of pancake lattices in the neighboring layers. A 
hexagonal pancake lattice in tlie nth layer (closed circles) is 
shifted with respect to the lattice in the (n— l)th layer (open 
circles). Vector s, eq.(34), is shown in a box. 

as a sum of smooth and oscillating parts: 

V|| = v|| -F(5v||, (35) 



FIG. 3: To the definition of vector g. When two pancake 
lattices in two neighboring layers are shifted with respect to 
each other pancakes got connected by Joscphson vortices in 
such a way as to minimize the length of the latter. The op- 
timal choice is shown by solid lines. Two other possible con- 
nections (broken line and dotted line) require longer vortices 
which, obviously, costs more energy. 
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It is non-zero only when is the reciprocal lattice vector. 
If qii belong to the reciprocal lattice, one prove 



giq||g(s„) _ giq||S„ 



(40) 



Thus, we have for the oscillating part: 



5v||(r,n)= (41) 



Po sr^ g(s«) 



^giq||U„(r) _ giq||U„_i(r)^ g- 



■iq||r 



where the sum runs over the reciprocal lattice vectors. 
The above expression is explicitly periodic with respect 
to both u„ and u„_i. This means that if one of the 
layers is shifted by a lattice vector the Josephson vorticity 
remains unaffected by such a shift. 

To calculate the energy we need the transverse compo- 
nent of the Josephson vorticity: 



(5v||T(r,n) 



Po qiuqii 



(42) 



„iq||U„(r) 



||U„_i(r)^ 



(Note: in 2D vector product is a scalar, not a vector.) 
As for V||, it is transverse as long as u„'s are transverse. 

We are in position now to write down the Josephson 
energy: 



Stt 



Hj = Hn + Hj2 + Hj3,{^3) 
Y^f /Cj(r-r',n-m)g(s„)g(s;,),(44) 



H 



32 



Stt ^ 



Hj3 



E ^j(qihO) 



[q|| X g(Sn)] 



9||(q||g(s„)) 

(2 - 2cosq||S„) , 
*°'°E/EMq|hn-m)x 
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[qii X sti^-JlIqii X g(«».)] 

?^(q||g(sn))(q||g(sm)) 



,«q||Un 



-iq||U„ 



-iqiiUm- 



where we used the notation 

= u„(r') - u„_i 



(r')- 



(47) 



The term Hji describes the kinetic energy of supercur- 
rents induced by the Josephson vorticity. The kernel 
/Cj(Ar, n) is non-zero for |Ar| <C Ac (see Appendix C) 
since such currents spread over large regions. Therefore, 
two Josephson vortices apart from each other interact 
via these currents. The second term is purely local. It 
describes the increase of the free energy of an interlayer 
Josephson junction due to insertion of vortices. The third 
term describes similar effect: the modification of the free 



energy of an interlayer Josephson junction by a vorticity 
in another layer. As we will show below, this effect is 
small. 

At low temperature, where it is possible to write 
g(s) « s, Hamiltonian Hj is responsible for the Joseph- 
son contribution to the tilt modulus C44 of the vortex 
lattice [11]. The first term corresponds to the B"^ piece 
of C44 while the second and the third terms contribute to 
the so-called single- vortex part (cx B) of the tilt modulus. 

At sufficiently high temperature we cannot approxi- 
mate g by s. Instead we must follow different type of 
analysis. The term 



/Cj(r-r',n-m)g(s„)g(s:„) (48) 



J r,r' 

may be disregarded right away if |n — to| > 2 since at 
low Josephson coupling we could neglect correlations of 
the displacement fields in different layers: (u„Uto) ~ 0. 

Thus, for \n — m| > 2 

g(s„)g(s;„) « g(s„)(g(s;„)) + g(s;,)(g(s„)) = o (49) 

since (g(s)) = 0. The latter statement is obviously true 
for g is an odd function. 
We have: 

iJji cx E / /Cjg(s„)g(s;) ^50) 

m,n 

-EE / {^Jgq„gp„e*''"^""-"''-^'-'''"^""-""-^^ 
- /Cjgq„gp||e'^ii(""+i-"")+'Pii«-<_i)| , 



PI, 7^0 



X (45) where notation (24) is used and gqn 



-g- 



are Fourier 



transform coefficients of the odd periodic function g(s). 
As one notices from this formula, Hji is a product of four 
exponents of the form exp(iq||u). It will be shown that 
Hj2 depends on two exponents only. At high enough tem- 
peratures Hji may be neglected in comparison with H]2- 
This is because (i?j2) is proportional to (exp(iq||U||))^ 
while (Hji) is proportional to (exp(iq||U||))* and the ex- 
pectation value of the exponent vanishes when T grows. 

One could adopt a more formal approach. As shown in 
Appendix B at sufficiently large T operator Hji becomes 
irrelevant in the renormalization group sense while Hj2 
remains relevant up until higher temperature. In this 
temperature interval it is possible to neglect H,n- its 
contribution to the thermodynamics is purely perturba- 
tive while the contribution of the relevant second term is 
singular. 

The term Hjs could be also disregarded for the follow- 
ing reason. The kernel 



/Cj(q) 



1 



Ac 



(51) 



(52) 
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may be simplified provided that the Josephson length is 
bigger than the lattice constant; 



Aj = 7dc > do- 
Namely, we can write: 

1 ll 



/Cj(q) 



1 



(53) 



(54) 



The term q^/^'^q^^ is smaller than unity as long as (53) 

is valid. Indeed, qz < tt/^c and g|| > a^^. This means 
that /Cj(q||,n) is smaller than /Cj(q||,0). Therefore, i7j3 
(corresponds to /Cj(q||, n)) is much less than ffj2 (corre- 
sponds to /Cj(q||,0)). 

Thus, we retain Hj2 only: 



Stt 



[q|| X g(Sn)] 



9ii(qiig(s„)) 



X (55) 



(2- 2cosq||S„) , 

= , \ 2 -(56) 

We rewrite this expression in a more compact way: 



3$2 



^ /'/i(u„(r)-u„_i(r)), (57) 



Ms) = E !l 



[qii X g(s) 



(qiig(s)) 



fl — cosqiis) . (58) 



Function h is dimensionless and periodic with the pe- 
riod of the hexagonal lattice. Its contour graph is shown 
on fig. 4. Due to its periodicity h can be expanded in a 
Fourier series: 



h{s) = /lo + E cosq||S. 



(59) 



Using this expansion it is possible to simplify the ex- 
pression for iJj at finite temperaturc;s with the help of 
the trick already discussed at the end of Sec. II: 



(60) 



3|/i,ol$^ 



2567r5A2 



^E_/^ E {l-cos(q||[u„-u„_i])}, 



|qiil=9o 



3lL 



-3.7, 



2567r5 



1.4 X 10"'', 



(61) 
(62) 



Here the sum over vectors q|| runs only over the shortest 

of them. The number hq^ is the Fourier coefficient of h{s) 
corresponding to such reciprocal vectors. This is the final 
equation for the Josephson energy of the vortex system. 



IV. COUPLING CONSTANTS 

Let us estimate coupling constants of our Hamiltonian. 
That way we can understand what are the largest interac- 
tion in our system and when the perturbative treatment 
of interactions is permissible. 

First, we calculate the elastic constant Sei- To find it 
we measure area in units of Pq ^. Thus: 



Eel 



647r2A 



2 ■ 



For BSCCO £ei = 182 K assuming: 

dc = 1.5 nm, 
Xab = 200 nm. 



(63) 



(64) 
(65) 



The Josephson constant corresponding to Hj2'- 

3$2 

-^= 256.5Akpo '"-'- ^''^ 
The coupling constant corresponding to Hji : 



(67) 



Using the estimate for the Fourier coeflacients rj l/2g'o 
we find: 



3$gao 



10247r4Aa6Ac 



(68) 



Ratios of these two coupling constants to the elastic en- 
ergy scale: 



gol _ 12|/igp| gp ^ 2 



(69) 
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FIG. 4: Contour plot of the function h{s) = h{x,y), defined 
by eq. (58) . Function h has a minimum at the origin and six 
maximums at the corners of the primitive cell. 
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and 



Eel 



6 an „ On 
IDTT^ Aj Aj 



(70) 



are both smaller than unity when (53) holds. This ex- 
plains the physical significance of (53). The latter is es- 
sentially a criterion for our vortex system to be viewed as 
a quasi-2D pancake gas. If (53) is violated the system is 
better described in terms of the Abrikosov vortices rather 
than pancakes. 

The interlayer magnetic interaction constant could be 
estimated directly from (29): 



3$§4 



-mag 



(71) 



The quantity 1/qodc is an estimate for number of layers 
coupled by the magnetic interaction with a given layer. 
The ratio of the magnetic and the Josephson constants 
is: 



-mag 



X2 



X2 



(72) 



The last result is important for it allows one to judge 
when the Josephson coupling could be neglected and 
when it must be retained. As one can see from the above 
formula, the ratio is independent of the magnetic field 
and, in that sense, it is a "material constant" . If we take 
for BSCCO: 



7 = 300, 



(73) 



then the ratio equals to 1.5. That is, even for such an 
extremely anisotropic material the Josephson coupling is 
of the same order as magnetic. 



V. UPPER BOUND FOR THE MELTING 
TEMPERATURE 

In this section we find the upper bound T„ on a melting 

temperature T^. The gist of the following calculations is 
that at sufficiently high temperature the interlayer cou- 
pling becomes irrelevant and the system could be thought 
of as a collection of 2D layers decoupled from each other. 
In these layers no long-range crystalline order is possible. 
The temperature found in this fashion is not necessary a 
true melting point for melting could occur at even lower 
temperatures through a first order phase transition. 

In certain situation our upper bound could serve as a 
reasonable estimate for the melting temperature. This 
happens if the Josephson coupling is sufficiently strong 
and the difference between the solid and the liquid phases 
of the pancake matter is small. When this difference is 
exceeded by the bare Josephson energy scale the phase 
transition is controlled by the Josephson energy. Since 
the latter vanishes at T„ the melting line lies close to 



We assume that the temperature is large so that 
Josephson term is irrelevant (see Appendix B). Only 
Hj2 part of the Josephson coupling and the magnetic 
coupling, eq.(32), have to be accounted for. In such a 
regime it is possible to get using simple perturbative 
argument which we are ready to present. 

The crystalline phase is characterized by a non-zero 
order parameter: 



A = (e"iii"") 



(74) 



With the help of A the interlayer interactions (32) and 
(45) could be written in a mean-field manner as 



cosq||(u„ - VLm) ~ Acosq||U„. 



(75) 



The above transformations reduces our multilayer prob- 
lem to a problem for a single layer with the Hamiltonian: 

-ffmf = ^ / ^|Vu„|^ - aA ^ cosq||U„, (76) 

" lq||l=9o 

a = Po(£j + Emag)- (77) 



It is presumed that a is small: 



(78) 



This guarantee the applicability of approximation (75). 

It is equivalent to (53). 

Let us transform the interlayer interaction term as fol- 
lows: 



- aAcosq||U„ « -aA^ ^1 - ^(qyu^)" 



aA- 



^(qiiu„)' = 2«"' 

a = a9^(cosq||U„)2 = a^^exp (-((qnu^J^^) . 



(79) 



(80) 



After such substitution the effective Hamiltonian be- 
comes quadratic. Thus, the problem could be solved ex- 
actly. However, when solving it we must keep in mind 
the self-consistency condition: 



((qii"J') 



T(q||k)2 



STrCee 



In^ 



P(C66fc2-f 5) 

In — z h const. 



((q||U„)^) + const. 



(81) 



The last equation has positive solution for ((qyu^)^) only 
when 



T<T„ = 



SttC, 
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V3$|4 
647r3A2,- 



(82) 



Above r„ the order parameter A is zero. 
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Since Xat is a function of temperature itself, to find T„ 
it is in fact necessary to solve the equation: 



T = 



647r3A2,(T„)- 



(83) 



Assume that the temperature dependence of Xab is given 
by a phenomenological formula: 



x-jiT) = xjmi-TyT^), 



(84) 



where the superconducting transition temperature Tg 
equals to 100 K for BSCCO. Eq.(83), thus, gives: 



T„,«62K. 



(85) 



Let us reiterate, that this result is "high-field" . Namely, 
it is derived under the assumption that (53) is satisfied. 
The latter inequality is equivalent to: 



„ ^0 



(86) 



Thus, there is no contradiction if at low fields {B < 
$o/Aj) the melting temperature is higher than T„ found 
above. 

We can also find the contribution of the interlayer 
(Josephson and magnetic) coupling to the total free en- 
ergy density of the pancake lattice: 

Finter/S = -a(cos(q|| u))^ = (87) 

f a \ 8nCee-<i'^T 



-aexp(-((qi|u)2)) = -Q 



66 



where the following expression for ((q||u)^) was used: 



((qyu)^) 



STrCee-q^T \ a 



The latter expression is a, trivial consequence of (81). 

Note that a/Cee ~ (ej + £niag)/£ei and qlT/SwCee = 
T/Tu- Thus, the formula for Fi„ter could be rewritten: 



Finter/S ~ — C( 



ma 



66 



1/(1-T/T„) 

2 „ I (89) 



Interlayer contribution to the free energy, as expected, 
vanishes at T = T„. This result is useful only if the 
lattice remains stable close to T„, that is, if is close 
to T,,. 



VI. DISCUSSION 

In this paper we proposed a description of the quasi- 
2D pancake lattice. The description is a version of an 
elasticity theory: the Hamiltonian is a functional of the 
displacement field only. Such a transparent structure al- 
lows easy generalization. The most obvious extension is 



for the case of non-zero pinning. This is an important 
direction for future studies. 

Another advantage of the method is that at hight mag- 
netic field (eq.(86)) all non- harmonic terms are small. 
Therefore, consistent perturbative analytical tools are 
applicable. 

At present we do not include dislocations and vacan- 
cies into our model. Yet, to certain extent the effect of 
both dislocations and vacancies may be accounted for 
even at this stage. Unlike the interlayer coupling, these 
two disorder the lattice. The interlayer coupling tries 
to eliminate them by binding into a topologically neu- 
tral pairs. In turn, dislocations and vacancies attempt to 
destroy the interlayer coupling by breaking translational 
order. Which side wins this competition could be de- 
termined through comparison of corresponding singular 
corrections to the free energy. At least conceptually, this 
is a straightforward task. 

The main goal of this paper is the derivation of the 
model. However, we'd like to discuss the most elementary 
of its consequences here. 

First of all, our model allows to compare quantitatively 
the magnitude of the magnetic and Josephson pancake 
interactions. We have seen above that even for BSCCO 
(7 ~ 300) both of these are of the same order. Therefore, 
models which neglect the Josephson interaction are only 
qualitatively accurate, at best. 

Next, let us briefiy discuss the application of formula 
(83) for pancake lattice melting in real layered super- 
conductors. Probably, it has no physical significance for 
BSCCO due to its extremely small Josephson coupling. 
The total interlayer coupling is also small since the mag- 
netic interaction is of the same order as the Josephson. 
The transition into the liquid phase occurs because the 
free energy of the liquid becomes less than the free energy 
of the solid. Small contribution of the interlayer interac- 
tion could change the melting line slightly but otherwise 
is insignificant. 

We would like to conjecture cautiously that (83) may 
be of relevance for Pb-doped BSCCO. This superconduc- 
tor has much smaller anisotropy 7 70 [12]. Presumably, 
this indicates that its Josephson coupling parameter is 
higher than that of pure BSCCO. 

In [13] Pb-doped BSCCO was studied by means of 
magneto-optic technique. A depinning transition was ob- 
served at T* « 54 K. The authors of [13] interpreted 
the transition as a pancake lattice melting. Remarkably, 
they see that this transition is field independent. This 
is a strong argument in favor of (83) being applicable 
since T„ does not depend on the magnetic field. Another 
field-independent melting transition, 2D dislocation un- 
binding, is an unlikely candidate for it occurs at much 
smaller temperatures. 

Pb-doped BSCCO has Tc « 91 K and 4 ~ 1.9 nm. In 
[14] Xab was reported to be of the order of 180 nm. Thus: 



Tu w 68K. 



(90) 



This is close to the transition at T* « 54 K reported in 
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[13]. The fact, that T,, and T* differs by about 25% is not 
by itsch discouraging since penetration depth \ab shows 
strong sample-to-sample variations [15]. It is very weh 
possible that they are responsible for the deviation of the 
experimental transition temperature from its theoretical 
estimate. Yet, the final judgment about the applicability 
of the presented theory is to be postponed until further 
investigation. 

In conclusion, we proposed a simple model for the pan- 
cake lattice. It could be used to study pinning and ther- 
modynamics of the pancake matter at sufficiently high 
fields and temperatures. 
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APPENDIX A: THERMALLY INDUCED 
JOSEPHSON VORTICITY 

In this Appendix we will try to estimate the charac- 
teristic size £,ioo-p(T) of the thermally induced Josephson 
vortex loops at given temperature T.We will also discuss 
what happens when ^loop grows bigger than cq. 

To achieve the first goal we calculate the energy of 
a single Josephson vortex loop -Bioop af' a function of 
its radius R and compare this energy to the temperature. 
We will assume that ^loop is determined by the condition 



-^loop (^loop ) — 



(Al) 



To proceed with the calculation we first find the Fourier 
transformation of the vorticity field of a single loop whose 
radius is R: 

V||(q) = / e^iii'"[z X dv\ = 2nRJi{q^\R) ^ ^"l (A2) 

Jr=R QW 



The energy is given by: 

^2 r 47r2ji2((7||i?) 



-E'loop(-R) = -^^^ I T 

All ,9. 1 



(A3) 



The integral in this formula could be transformed as fol- 
lows: 



/ T 



AT^^JliqnR) 



-m 



2 / ' arctg 



7rA„ 



q\\dq\\, 



where ^ is the superconducting coherence length. After 
a substitution x = gyi? the integral becomes: 

j Jl{x)Bxctg y-y-^ j dx. (A5) 



Xab^cR 

In deriving the above equation we used the approxima- 



tion ^1 -\- AcQy ~ Ac^ii . This is accurate for gy ^ ""^ 

or, equivalently, R^Xc- 

Finally, confining ourselves to the range ^ <C i? <C Aj 
we write: 

2 fTTR\ f°° .dx TT 



Therefore: 



El 



f rV 87r3„ f R 



(A7) 



where T„ is given by (83). Using (Al) one determines: 

(A8) 

At T = T„ which is the highest temperature at which the 
model proposed in this paper is still applicable we have 



600P ^ 0.08Aj. 



(A9) 



This give us the upper bound for the characteristic size 
of the thermally induced Josephson vorticity loops. 

As it was discussed in the body of the paper we would 
like to have ^loop o-o (eq.(33)). Consequently, 



O.OSAj < ao. 



(AlO) 



'='^11, 



It is not difficult to notice that the latter condition is not 
always compatible with the requirement ao ^ Aj which 
has to be imposed to observe the quasi-2D effects of the 
pancake physics (eq.(53)). What does this mean for the 
analysis undertaken in this paper? 

The answer to this question is that it is possible to 
recover virtually all of the physics described in the paper 
even if (33) is not satisfied as long as ^loop ^ ^o- (The 
quantity ^0 was defined in the paper as a size within 
which the 2D pancake lattice could be viewed as almost 
ideal.) Let us examine how this is done. 

Presumably, the violation of (33) might lead to prob- 
lems with (57). It was calculated under assumption that 
it is possible to define a single-vortex line tension on a 
scale ~ ttQ. When vorticity fluctuation scale ^oop exceeds 
ao it becomes difficult to assign a well-defined line ten- 
sion to a Josephson vortex. Thus, the analysis at the 
scale of ao fails. 

It is, however, possible to study the system at the scale 
of ^0- At this scale we need not to be concerned with the 
properties of a single Josephson vortex. Instead, we will 
think in terms of average local vorticity which pierces 
interlayer Josephson junction. 
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Imagine first the situation when pancake lattices in 
the nth and (n — l)th layers coincide at r = ro. This 
implies that for any r, |r — ro| < ^o, the average Joseph- 
son vorticity between these two layers is absent. Thus, 
the local density of the Josephson energy is zero. Now 
we shift the lattices with respect to each other: s„ = 
u„(r) — u„_i(r) ^ 0. The average vorticity (v||(s)) is 
no longer zero as well. We have to disentangle its contri- 
bution from contribution of the fluctuation-induced vor- 
ticity 5v|| in the situation where fluctuations locally are 
strong: (|(5v|||) > |(v||)|. To do so we coarse-grain the sys- 
tem description up to the scale ^loop- After this coarse- 
graining (5v|| K. and the bare value of the penetration 
depth has to be replaced by Ac- (The latter is the 
penetration depth measured experimentally.) Thus, at 
the scale ^ Cioop the vorticity fluctuations could be 
neglected and it is possible to work with the effective 
Josephson coupling instead of including fluctuations ex- 
plicitly. 

The "mean-fleld" vorticity (v|| (r)) is a periodic func- 
tion of s„. This implies that the pancake lattice energy 
density due to Josephson coupling is a periodic function 
of s„ as well: 



it docs not rely on ^loop ^ condition. The advantage 
of the latter, however, it its ability to give the numerical 
constant accurately. The fact that numerical coefficients 
of (45) and (A14) are very close is pure luck. Obviously, 
the estimates like (A13) have "order of magnitude" pre- 
cision only. Numerical constant in (A14) should not be 
taken too seriously. 

We re-derived Hj^ above. Operator Hj\ may be re- 
derived in the similar fashion. 



APPENDIX B: SCALING DIMENSIONS 

In this Appendix we determine the scaling dimension 
(SD) of the operators (50) and (60). We start by solving 
an auxiliary problem: we calculate SD of an exponent 
exp(iq||u). The easiest way to find SD is to evaluate the 
correlation function: 



^e'Hi"We~''iii = e~5((qiil"-"'])''). (Bl) 

The correlation function of the displacement field is 
found: 



-i((q||[u-u']f> 



^J^E/^J/^o E [l-cos(q||S„)]. (All) ((q|,[u-u'])2)=|^^g^(l-e^'^(--')). (B2) 



|qiil=9o 

This expression is completely general: it is the lowest qy 

terms of periodic function Fourier scries. 

The final question now is the estimation of £j. The 
latter, however, is fairly trivial. As it was discussed in 
the body of the paper iJj2 comes from the increase of 
the interlayer Josephson junction energy due to non-zero 
(v||(r)). Thus, we have match (All) and the interlayer 
Josephson junction free energy: 



The integral over k was evaluated in [16] (eq. (9.3.33) 
and problem (9.9) of the latter reference): 

l^'^il- e'"^-) = l-S,^ In (|r|/ao) + . . . , (B3) 

where ellipsis stand for terms which remain bound at 
|r| — > oo. Consequently: 



^j(ro) 



$2 



167r3A2(i, 



[1 - cos{</>„(ro) - (/)n_i(ro)}](A12) 



((q||[u-u'])^) 



47rC66 



ln|r-r'|/ao + ---- (B4) 



When s„ = both are zero. When |s„| ~ Oo/2 the 
Josephson junction energy J^j is no longer zero but rather 
equals to <I>Q/167r''A^(ic. It is because the cos term in the 
above expression effectively vanishes due to fast oscilla- 
tion of the phase difference. This implies that 



Therefore: 



2T/T„ 



^giq||u(r)g-iq||u(r')^ ^ 



(B5) 



From this equation we can read off SD of the exponent: 



sjpo max < ^ [1 - cos(q||s) 

. |qiil=9o 



..(A13) 



(B6) 



Since maximum of the sum in this formula is equal to 9 
the coupling constant ej is: 



£jPo 



1447r3A2dc 



2.2 X 10" 



(A14) 



Equations (All) and (A14) reproduce (45) up to a nu- 
merical constant of order unity. 

The derivation of (45) developed in this Appendix is 
more general than that found in the body of the paper for 



where T„ is given by (82). Thus, SD of (60) is equal to 
2T/T„ since it is a product of two exponents. This SD 
becomes equal to 2 when T = T„. At T > opera- 
tor (60) is irrelevant, otherwise it is relevant. As usual, 
this mark the transition point from ordered to disordered 
phase driven by (60). Eq.(82) establishes the same fact 
derived by other means. 

As for (50), its operator density is proportional to: 

/ ^g(sW)g(s(r + R||)). (B7) 
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Here the factor ^/R\\ comes from the kernel /Cj. The 
kernel /Cj is calculated in Appendix C. 

SD of g is 2T/T„ as it is obvious from (50), SD of the 
integration measure is -2, SD of 1/-R|| is 1. Totally, we 
have 



= 4T/T„ - 1. 



(B8) 



Equating this SD with 2 we discover that for T > 3T„/4 
operator (50) is irrelevant. Therefore, is is permissible 
to neglect it for these temperatures. At the same time, 
operator (60) remains relevant up until the temperature 
has grown to T„. Thus, the discussion of the Sect.V is 
valid for 



3T„ 



<T <Tu. 



(B9) 



APPENDIX C: KERNEL /Cj 



In this Appendix we calculate the kernel /Cj(r, 0) re- 
quired for eq.(50). 

We must evaluate the integral: 



/Cj(r,0) = 



k„y_„/rf^ 277 1 + A2fc2+A2,fc2 



(CI) 



arctg 



In most layered superconductors the screening length Ac 
is very big. Therefore, wo are interested in distances 
r <C Ac. For such r it is possible to neglect unity under 
the square root and write: 



J- / 

^a6Ac Jkii 



-arctg 



(Ajfcii) 



Next we perform the angular integration: 
1 



(C2) 



l/oo 



Jo(a^)arctg | — | ax. 



27r2A, 



1 r 

^abKr Jo 



A J a; 



The upper limit in this integral is assumed to be large: 
r > ao. 

If r <C Aj then arctg(7rr/Ajx) could be approximated 
by nr/Xjx for almost whole integration interval. The 

contribution of the interval where this approximation is 
invalid could be evaluated separately: 



/Cj 



27r2A 



1 f r''/^-' / Trr\ 



+ 



/ Mx)^ — dx > . 



The first integral is 0(7rr/Aj). We neglect it. The second 
integral must be integrated by parts: 



/ -^^dx w - In — + / Ji{x)\nxdx. (C5) 

J-rrr/X, ^ Aj Jq 



f rr/ao 
7rr/Aj 

Consequently 



/Cj 



1 



, 7rr 

In h const. 

Aj 



for r «: Aj. (C6) 



27rA24 

In the opposite limit r ^ Aj we define a function: 



f{x) = [ dyJo{y), 
Jo 



(C7) 



and write: 
/Cj« 



1 



2Tr^ Xab\ 



c' Jo 



r/ao 



arctg [ — ]df= (C8) 



27r2Aa6Acr VAj 



7rr 



f{x)dx 



a;2 + (7rr/Aj)2- 



(C9) 



The integral may be transformed in the following man- 
ner: 



f(x)dx 



f{x)dx ^ 
Jo x^ + (7rr/Aj)2 io 

rr/ao 

Jl X"^ + '-r-r/\.\2 

Jo x2 + (WAj)2+^("')A ^ 



+ (7rr/Aj)2 ^ x^ + inr/Xj^ 

f{x)dx 
x^ + (7rr/Aj)2 

dx 



-(CIO) 



2 + (7rr/Aj)2- 



In the last line the first integral is smaller than the second 
one. Indeed, the first one is 0((Aj/r)2) while the second 
is ©(Aj/r). Therefore: 



rr/ao 

Jo 



f{x)dx 



/r/ao 



a;2 + (7rr/Aj)2 

dx Aj 
a:2 + (7rr/Aj)2 ~-'^°°^2r" 



(Cll) 



Note that: 



/(oo) = / Jo{x)dx = 1. (C12) 

JO 



Thus, for r <C Aj we have: 



/Cj 



1 



47rAabAcr 



(C13) 
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